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ABSTRACT 


Drag  measurements  on  spheres  in  high  speed  transition  regime  flow  are 
presented.  The  spheres  are  suspended  e I ectromagneti ca I ly  in  the  low  density 
flow  field  from  a  jet  expanding  freely  from  a  small  sonic  nozzle  into  a 
vacuum.  This  arrangement  provides  sting-free  measurements  under  hypersonic 
conditions.  The  current  in  the  control  coil  of  the  electromagnetic  balance 
is  proportional  to  the  applied  force  and  provides  a  sensitive  determination 
of  the  small  forces  encountered. 

Data  were  taken  using  nitrogen  and  argon  gases.  A  number  of  nozzle 
and  sphere  sizes  were  employed  covering  a  Knudsen  number  range  of  0.05  to 
5.  The  results  exhibit  a  smooth  increase  in  +he  transition  regime  drag 
coefficient  toward  the  free  molecular  limit  for  diffuse  reflection  and 
complete  thermal  accommodation. 

Comparison  is  made  with  the  available  experimental  results  ov  other 
techniques.  Improved  repeatability  and  an  extension  of  range  of  flow 
parameters  is  obtained  with  the  present  methods.  The  data  are  compared 
with  current  near-tree  molecular  flow  theories  and  the  modified  Krook 
solution  of  Willis  is  found  to  give  the  best  agreement  with  the  experi¬ 
mental  results. 
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SECTION  i 
INTRODUCTION 


The  primary  objective  of  this  study  is  to  obtain  sphere  drag  data 
of  improved  precision  in  transition  flow  at  high  speeds  with  emphasis  on 
near-f ree-mol ecu le  flow  conditions.  The  results  presented  in  this  report 
were  obtained  using  free  magnetic  suspension  of  the  model  in  a  free  jet 
flow  field.  These  techniques,  while  introducing  some  experimental  in¬ 
conveniences,  eliminate  two  basic  objections  to  wind  tunnel  testing  on 
three  dimensional  models  at  very  low  density;  viz.  the  uncertain  sting 
or  mechanical  support  effects,  and  excessive  pumping  requirements.  In 
addition,  the  unique  suspension  system  has  a  sensitivity  which  provides 
a  degree  of  precision  and  reproducibility  which  has  not  heretofore  been 
available  for  the  measurement  of  the  small  forces  encountered  in  low 
density  flows. 

The  transition  regime  of  gas  dynamics,  between  continuum  and  free 
molecular  flows,  has  been  the  object  of  much  theoretical  and  experimental 
investigation  in  recent  years.  Theoretical  analysis  of  transition  flow 
is  complicated  by  the  presence  of  the  collision  term  of  the  Boltzmann 
equation.  When  the  flow  is  truly  free  molecular  some  aerodynamic  problems 
can  be  solved  since  only  molecule-surface  collisions  s^e  present.  Even 
in  this  limiting  case  difficulties  arise  since  +here  arc  unresolved 
questions  concerning  gas-surface  interactions.  In  transition  flow, 
intermolecu lar  collisions,  in  addition  to  mo  I  ecu! e-surface  collisions, 
complicate  fhe  problem  beyond  mathematical  comprehension  unless  simplifying 
assumptions  are  made.  To  date  >t  has  been  necessary  to  impose  severe 
restrictions  on  the  analysis  which  delegates  much  responsibility  to  ex¬ 
perimental  results.  Gross  properties  such  as  drag  provide  one  Input  to 
the  solution  of  transition  flow  problems  and  sphere  drag  measurements 
provide  simple  three  dimensional  aerr.  ynarnic  data.  The  geometry  is  con¬ 
venient,  yet  it  provides  a  more  re?  'stic  model  for  analysis  than  two 
dimensional  shapes. 


Satisfactory  experimental  data  are  available  for  two  dimensional 
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geometries  in  transition  flow.  ’  Problems  of  size  and  mounting  in  this 
case  are  such  as  to  be  readily  adaptable  to  more  conventional  experi¬ 
mental  techniques.  However  mechanical  mounting  of  3-d i mensi ona I  body  such 
as  a  sphere  is  difficult.  If  the  body  is  made  sufficiently  large  to 
reduce  sting  effects,  it  is  difficult  to  achieve  high  Knudsen  Nos  over 
the  range  of  high  velocities  which  is  of  qeneral  practical  interest. 

On  the  other  hand,  reducing  body  size  leads  to  difficulty  in  meansuring 
the  small  forces  which  exist  in  low  density  flow. 

Electromagnetic  model  suspension  eliminates  an/  sting  or  support 
effect  and  allows  the  accurate  determination  of  small  forces.  The  use 
of  this  technique  has  resulted  in  sphere  drag  measurements  with  little 
scatter  and  good  reproducibility.  The  minimum  sphere  size  that  can  be 
used  is  limited  only  by  mechanical  handling  problems  and  the  optical 
technique  used  with  the  support  servo  system.  With  these  imitations 
high  speed  transition  flow  relative  to  the  sphere  can  be  produced. 

The  technique  used  to  provide  the  low  density  flow  field  employs 
the  freely  expanding  jet  from  a  converging  nozzle.  These  jets  have  beer, 
the  object  of  much  investigation^  and  many  of  their  properties  are  known. 

A  high  velocity  low  density  flow  is  produced  when  the  gas  expands  into 
a  vacuum  and  the  core  flow  is  isolated  from  the  surroundings  when  the 
proper  operating  conditions  are  observed.  Properties  on  the  jet  axis 
can  be  calculated  with  considerable  confidence  and  a  reasonable  range 
of  density  and  Mach  number  is  available.  The  problems  encountered  due 
to  the  expanding  nature  of  the  flow  are  outweighed  by  the  difficulty  and 
expense  of  obtaining  similar  flow  conditions  by  other  means. 

Using  a  combination  of  the  magnetic  suspension  and  free  jet  tech¬ 
niques,  sphere  drag  measurements  have  been  made  in  the  range  from  transi¬ 
tion  flow  to  nearly  free  molecular  conditions  in  Argon  and  Nitrogen 
(.01  Km  <  5).  The  theoretical  inviscid  Mach  Nos.  ranged  from  8.7  to 
17.8  for  nitrogen  and  from  15.0  to  25.6  for  Argon.  The  data  are  compared 

X 

See  for  example  references  4  and  5. 
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with  the  available  experimental  results  of  other  techniques 
as  a  suggested  semi -emp i r ica I  fit  and  theoretical  models. 
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SECTION  II 

THE  WIND  TUNNEL  BALANCE  SYSTEM 


A.  The  Magnetic  Suspension  System 

Free  electromagnetic  suspensions  controllable  in  one  dimension  have 

been  in  use  since  the  1 930 ’ s .  However,  systems  controllable  in  three 

dimensions  having  a  wide  applicability  to  problems  of  aerodynamic  interest 

are  relatively  recent.  In  ! 957  the  concept  of  a  three  axis  suspension 

system  in  which  the  forces  on  the  supported  object  were  automatically 

resolved  into  three  Independent,  mutual iy-perpendicular  components  was 

devised  by  Parker.^  Experimental  verification  of  the  concept  was  reported 
7  8 

in  1959.  ’  Development  in  a  form  suitable  for  use  as  wind  tunnel  balance 

9 

was  completed  by  1964  and  the  present  design  became  operable  in  1965. 

This  design  as  well  as  fhe  genjral  theory  of  operation  has  been  described 

■  .  ,*  4.  10,11,12“ 

in  the  htera+ure.  '  ' 

The  actual  coil  system  used  in  the  present  apparatus  does  not  follow 
any  ideal  configuration  precisely,  but  rather  utilizes  a  compromise  ar¬ 
rangement  conceived  for  the  specific  application.  The  major  criterion  was 
maximum  force  capability  in  the  vertical  direction  while  maintaining  a 
reasonable  tunnel  size.  The  system  provides  independent  three  dimensional 
control  and  results  in  a  linear  re  I ai ionshi p  between  applied  force  and 
coi !  current. 

I .  Coi I  Sets 

The  coil  system  is  shown  in  Figure  2.1.  A  pair  of  large  coils,  A,  is 
used  +0  provide  a  magnetic  moment  per  unit  volume  in  the  model.  Gradient 
coil  pairs,  B.  C,  are  located  within  the  core  of  the  main  field  coils  and 
produce  forces  on  the  model  in  accordance  with  the  general  relation 

A. 

dF  =  (dm  •  V)B  (2.1) 

An  analysis  of  these  magnetic  forces  as  a  function  of  gradient  coil 
geometry  is  given  by  Jenkins  and  Parker,  and  there  is  a  wide  variety  of 
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FIGURE  2.1  Coil  Configuration 


geometrical  arrangements,  including  the  one  used  here,  which  wi 1  I  produce 
forces  readily  resoluable  into  three  mutually  orthogonal  components. 

The  airflow  is  directed  vertically  downward  and  hence  the  current  on 
the  vertical  gradient  coils  is  directly  proportional  to  the  drag  on  the 
model.  In  the  case  of  the  sphere,  the  gradient  coil  sets  producing  forces 
in  the  horizonal  plane  are  used  only  to  provide  lateral  stability. 

A.  Servo  Control 

The  inherent  instability  of  a  magnetic  suspension  system  requires  a 
servo  control  system  for  each  of  the  three  degrees  of  freedom.  The  "error" 
signal  to  actuate  this  control  is  providea  by  two  light  beams  at  90°  in 
the  horizontal  plane.  After  passing  over  the  model,  each  beam  falls  on 
a  separate  pair  of  photodiodes  located  symmetrically  about  the  axis  of 
the  sensor  system.  Thus  any  motion  of  the  sphere  from  its  null  position 
will  generate  changes  in  output  of  the  photo  diodes.  When  properly 
compared  these  outputs  will  yield  appropriate  signals  in  each  of  the 
orthogonal  directions.  With  proper  processing  these  signals  are  used  to 
vary  the  currents  in  the  respective  gradient  coil  sets  so  as  to  maintain 
the  sphere  at  the  null  position.  A  measure  of  the  current  change  required 
is  a  measure  of  the  force  causing  the  displacement. 

3.  Force  Cal i brat  ion 

A  calibration  was  made  to  verify  linearity  and  account  for  any 
gradient  that  might  occur  due  to  coil  asymmetry.  Initially  the  giadient 
coil  pair  symmetry  plane  was  located  as  accurately  as  possible  with  a 
differential  gaussmeter.  Then  one  of  the  main  field  coils  was  shunted 
until  the  symmetry  plane  of  these  coils  coincided  with  the  zero  field 
plane  of  the  verfical  gradient  coils.  This  fixed  the  position  of  sphere 
support  in  the  field. 

For  a  calibration  of  force  against  gradient  coil  current  nonmagnetic 
weight  must  be  added  to  the  spheres  so  that  the  magnetic  moment  induced 
by  the  main  field  at  fixed  field  current  remains  constant  as  the  lead  force 
is  increased.  In  addition,  any  such  dead  weight  must  be  added  symmetrically 
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to  avoid  applied  torque  since  the  spheres  have  a  preferred  orientation  in 
the  field.  Finally,. the  resulting  diameter  cannot  be  significantly  larger 
than  the  original  sphere  and  be  compatible  with  the  optical  sensing  system. 

The  techniques  used  to  produce  this  calibration  weight  included  the 
application  of  various  coatings  to  the  spheres  and  the  attachment  to  the 
spheres  of  small  brass  rods  in  line  with  the  magnetic  moment. 

A  number  of  calibration  runs  were  performed  over  the  period  that 
the  data  were  taken.  Both  I  mm  and  1/8  inch  diameter  spheres  were  used  and 
the  results  were  repeatable.  The  last  such  calibration  is  shown  in  Figure 
2.2.  Each  point  represents  a  number  of  "weighings"  with  the  same  calibra¬ 
tion  sphere.  Eleven  different  spheres  were  user  »nd  both  operational  modes 
of  zero  and  25  mv  ^ixed  gradient  were  checked.  The  data  were  fitted  to  a 
straight  line  with  a  standard  deviation  of  0.IC8  mv  and  a  probable  error  of 
±  0.0723  mv.  The  y  inte-cepf  was  -  0.65  and  the  sphere  drag  data  includes  this 
force  correction.  The  probable  error  is  ±  0.17$  at  the  sphere  weight  (no  load 
condition)  and  since  all  forces  were  measured  as  the  sum  of  drag  and  weight 
(within  the  actual  calibration  range)  the  probable  error  in  any  total  force 
measurement  would  be  <.0.17$.  The  repeatability  with  the  "clean"  spheres 
used  for  the  drag  measurements  was  i  0.2$  of  sphere  weight  and  the  noise 
level  was  generally  less  than  0.2$  under  operating  conditions.  The  calibration 
was  checked  by  "weighing"  each  sphere  under  zero  gas  load  before  each  run  and 
after  those  runs  during  which  the  sphere  was  not  lost.  It  should  be  noted 
that  the  calibration  includes  the  least  accurate  force  measurements  made  due 
to  the  difficulty  of  eliminating  irregularities  in  the  calibration  models. 

Hence  the  balance  is  believed  to  be  more  accurate  than  is  illustrated  by  the 
data. 

The  basic  models  were  chrome  steel  alloy  spheres  and  their  masses  and 
sizes  are  given  in  Table  I. 
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FIGURE  2.2  Balance  Calibration 


TABLE  I 


SPHERE  MASS  VALUES 


Sphere  Diameter 
i/8  In. 

I  mm 

0.025  in. 


Mass  (mg) 

130.0  ±  0.2 

4.8  ±  0.01 

1.045  ±  0.010 


B.  The  Wind  Tunnel 

The  free  jet  used  as  a  flow  field  expands  from  a  small  converging 
nozzle  into  an  eight-inch  diameter  vertical  lucite  cylinder.  This  section 
passes  through  the  center  of  the  magnet  arrangement  and  contains  optical 
windows  for  passing  servo  system  light  beams.  The  cylinder  is  evacuated 
by  a  Roots  615  booster  pump  in  series  with  a  Stokes  412  H  Mechanical  pump. 
Bellows  were  provided  to  reduce  the  transmission  of  vibration  to  the  sup¬ 
port  system. 

This  system  has  a  flow  capability  of  about  1200  cfm  of  dry  gas  and 
the  tunnel  pressure  was  about  10  J  torr  under  no  gas  load.  The  character¬ 
istics  of  the  system  are  such  that  the  background  pressures  rises  as  the 
source  pressure  is  increased.  The  ratio  of  source  to  background  pressure 
remains  nearly  constant  for  a  given  nozzle  under  test  conditions. 

It  is  necessary  to  provide  for  alignment  of  the  nozzle  with  the 
sphere  and  for  the  positioning  of  the  nozzle  relative  to  the  sphere, 
since  the  flow  properties  in  the  jet  are  dependent  upon  the  axial  distance 
from  the  nozzle  exit. 

Also,  a  semi-automatic  sphere- I oader  was  devised  so  that  models 
can  be  placed  in  support  without  opening  the  flow  system.  This  is  of 
considerab'e  practical  convenience. 
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C.  Experimental  Procedure 

Selected  spheres  were  cleaned  and  Inserted  into  the  loading  chamber. 

The  loading  device  was  placed  in  the  vacuum  chamber  and  the  system  sealed. 
The  system  could  be  evacuated  within  a  few  minutes  to  aboul  10  torr. 
Lateral  and  vertical  nozzle  adjustments  were  made  and  the  desired  nozzle- 
sphere  separation  was  established. 

With  the  sphere  in  position  the  main  current  was  increased  1o  allow 
the  sphere  to  rotate  to  its  preferred  orientation  in  the  field.  The 
gradient  coil  currents  were  activated  and  the  vertical  current  was  increased 
until  the  sphere  lifted  slightly  from  the  loading  device.  The  loader  could 
then  be  lowered  completely  out  of  the  test  section  region.  Spheres  would 
remain  suspended  in  this  mode  under  no  gas  load  for  long  periods  of  time 
without  noticeable  drift  or  further  adjustment.  Recording  the  sphere 
support  current  (in  the  vertical  gradient  coil)  under  no  gas  load  condi¬ 
tions  provided  a  force  calibration  check  for  each  tes-f  run. 

Flow  was  established  by  admitting  gas  to  the  nozzle  from  a  dry  gas 
bottle  through  a  two  stage  regulator  and  needle  valve.  A  more  sensitive 
nozzle  alignment  could  then  be  performed  by  adjusting  the  nozzle  position 
for  zero  lateral  force  on  the  sphere.  Local  properties  could  be  varied 
by  adjusting  the  nozzle-sphere  separation  or  by  varying  the  source  pressure. 

A  data  run  consisted  of  incrementally  increasing  the  source  pressure 
and  recording  the  support  current  at  each  increment  with  the  nczz 1 3-sphere 
separation  remaining  fixed. 

Ideally,  the  source  pressure  was  increased  uni  i  I  the  maximum  force 
capability  of  the  balance  was  approached  and  then  decreased  fo  zero  to  pro¬ 
vide  a  final  check  on  the  force  calibration.  (In  practice,  however,  the 
sphere  was  often  lost  during  the  run.)  The  model  position  in  the  field 
was  carefully  maintained  by  adjustments  during  the  run. 

Gas  could  be  admitted  directly  into  the  vacuum  chamber  to  vary  the 
pressure  level  at  constant  source  conditions.  The  sensitivity  of  the 
force  measurement  to  the  chamber  pressure  level  or  the  pressure  ratio 
was  checked  for  each  set  of  nozzle  conditions. 
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SECTION  III 


THE  FREE  JET  FLOW  FIELD 


A.  Jet  Geometry 

The  characteristic  dimensions  of  the  free  jet  Lave  been  determined 
by  Bier  ana  Schmidt4  using  Schlieren  photographs  and  by  various  authors^ 
using  other  techniques  and  with  a  number  of  gases.  These  dimensions  are 
a  function  of  the  nozzle  size  and  the  ratio  of  source  pressure  to  back¬ 
ground  pressure,  P  ,  outside  the  jet.  The  core  flow  is  contained  within 
oblique  shocks  and  the  nearly  normal  shock  or  Mach  disc  downstream.  When 

the  proper  ratio  of  stagnation  to  background  pressures,  P  /P  ,  is  maintained 

o  c 

and  the  background  pressure  is  sufficient  for  the  shocks  to  be  well  formed, 
the  core  flow  is  independent  of  the  chamber  pressure  and  may  be  described 
by  isentroplc  assumptions. 

The  most  important  dimension  when  using  the  jet  as  a  flow  field  is 
the  axial  distance  to  the  Mach  disc.  This  may  be  expressed  as5 


V 

d  ' 
n 


1/2 


indepenuent  of  the  gas  used  for  15  <_  PQ/PC  17,000.  The  absolute  dimensions 
of  the  jet  flow  field  are  relatively  fixed  for  a  given  pumping  plant.  For 
the  wind  tunnel  balance  system  the  Mach  disc  was  about  32  mm  or  more 
downstream  for  bo+h  nitrogen  and  argon. 


B.  Flow  Variables 

Since  the  flow  properties  in  the  core  flow  vary  along  the  jet  axis, 
the  local  conditions  at  the  model  can  be  altered  by  adjusting  the  nozzle- 
sphere  separation  as  well  as  the  source  pressure.  Variables  character¬ 
izing  the  flow  field  can  be  calculate  using  familiar  gas  dynamic  concepts. 
A  Knudsen  number  defined  as 

*|* 

See  bib  I iography. 


(3.2) 


is  easily  obtained  if  It  is  assumed  that  p>.  =  Constant  on  the  jet  axis. 
This  implies  that  the  gas  collision  cross-sect  ion  can  be  evaluated  at 
stagnation  temperature.  The  assump+ion  may  not  be  unrealistic  since 
the  sphere  is  approximately  at  stagnation  temperature.  (This  is  dis¬ 
cussed  in  Chapter  V).  The  Reynolds  number  behind  a  normal  shock, 


Res  = 


p2V2d 


(3.3) 


can  also  be  evaluated. 

Familiar  nondimenslonal  variables  based  on  local  free  stream  condi¬ 
tions  require  a  knowledge  of  the  local  temperature  and  viscosity.  Due  to 
the  nonequilibrium  nature  of  the  flow  at  the  large  distances  downstream 
encountered  here,  different  temperatures  may  be  obtained  perpendicular 

4* 

and  parallel  to  the  jet  axis.  This  is  further  complicated  by  the 
"freezing"  of  the  parallel  temperature  at  certain  values  of  PQdn  while 
the  perpendicular  temperature  continues  to  decrease  in  a  manner  similar 
to  the  isentropic  prediction.  Even  if  tho  manner  of  combining  these 
values  were  clear,  viscosity  data  at  the  resulting  low  temperatures  is 
open  to  question.  However,  certain  comparisons  require  a  free  stream 
Knudsen  number  and  for  these  purposes  an  expression  based  on  the  Leonard* 
Jones  potential  model  will  be  used.  The  resulting  Knudsen  number  is 


Kn  =  - 1 -  (3.4) 

/2  n  7tp  ,2°d 
00  d  s 

where  ft  is  tabulated-in  Hirschfelder,  Curtis  and  Bird16  and  is  the 
hard  sphere  molecular  diameter.  The  value  ft  is  obtained  by  extrapolation 
and  based  on  the  isentropic  temperature  value  for  convenience.  It  should 

*1* 

'See  for  example  Muntz15. 
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be  noted  that  in  addition  to  the  fact  that  such  a  calculation  is  at  best 
very  approximate,  it  is  not  at  all  clear  that  this  is  a  reasonable  parameter 
for  correlating  drag  data.  In  any  event,  calculation  of  ft  based  on  typical 
freezing  values  of  T| |  do  not  significantly  alter  the  results  used  here. 

C.  Effective  Nozzle  Size 

Nozzle  throat  viscous  effects  are  important  at  low  densities  as  has 
been  indicated  by  Ashkenas  and  Sherman,5  Reis,17  Lefkowitz,18  and  others. 

The  jet  flow  should  scale  as  some  effective  nozzle  diameter  which  is  a 
function  of  the  nozzle  (or  orifice)  shape  as  well  as  the  throat  Reynolds 
number.  Although  it  Is  not  certain  that  the  correct  effective  size  is 
that  determined  from  mass  flow  discharge  coefficients,  it  appears  that 
the  use  of  such  a  diameter  would  be  more  realistic  than  using  the  actual 
throat  size  at  the  low  throat-Reynol ds  Nos.  encountered  in  the  sphere 
drag  experiments. 

Thus  mass  flow  measurements,  taken  with  three  nozzles  (each  having 
a  different  profile  were  used  to  calculate  the  effective  nozzle  diameter 
as  a  function  of  Reynolds  number  assuming  sonic  conditions  at  the  throat. 

The  effective  diameters,  normalized  to  the  actual  throat  diameter,  are 
shown  in  Figures  3.1  and  3.2  for  nitrogen  and  Figure  3.3  for  argon. 
Variations  in  geometry  bring  about  oitferent  results  for  each  nozzle. 

The  1,1  mm  diameter  nozzle  had  a  smoothly  varying  profile  approximately 
exponential  in  shape.  The  2.1  mm  nozzle  was  conical,  fairing  iuio  a  straight 
section  at  the  throat  and  the  0.95  mm  nozzle  was  a  relatively  thick  orifice. 

This  effective  nozzle  size  was  used  to  calculate  the  jet  flow 
properties.  It  is  interesting  to  note  that  the  I . I  mm  nozzle  results  and 
the  orifice  results  are  in  agreement  with  the  values  obtained  by  Ashkenas5 
using  impact  pressure  measurements  (there  is  3  slight  discrepancy  in  the 
high  Reynolds  number  asymtotes  indicated).  Lefkowitz18  also  found  agreement 
in  effective  nozzle  size  implied  by  impact  pressure  measurements  and  mass 
f low  va I ues. 
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FIGURE  3.1  Mass  Flow  Measurements  with  Nitrogen  (shaped  nozzle) 
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FIGURE  3.2  Hass  Flow  Measurements  with  Nitrogen  (tube;  orifice) 
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FIGURE  3.3  Mass  Flow  Measurements  with  Argon  (shaped  nozzle) 


0.  Flow  Gradients 


The  gradients  resulting  from  the  expanding  nature  of  the  flow  must  be 
considered  if  a  correlation  is  to  be  made  with  results  obtained  in  a  uniform 
stream.  The  effect  of  these  gradients  for  drag  purposes  can  be  characterized 
by  the  variation  in  density  since  the  velocity  is  essentially  constant  at 
large  nozzle-sphere  separations. 

Ashkenas5  gives  an  expression  for  the  radial  variation  of  density  that 
predicts  the  characteristics  calculation  within  about  3%.  This  equation  is 


pTOT^T  =  cos*9cos2<|£.> 


(3.5) 


where  r  is  the  radial  distance  perpendicular  to  the  jet  axis  and 


6  =  tan 


-I 


x  -  x 


(3.6) 


The  constant  <j>  has  the  v^lue  of  1.662  for  nitrogen  and  1.365  for  argon  and 
x 

■—  is  0.40  and  0.075  respectively.  A  calculation  for  the  worse  case  en- 
n 

countered  in  these  measurements  (largest  sphere,  smallest  s/d^,  argon) 
shows  that 


p (r,x) 
p (o,x) 


0.96 


(3.7) 


at  the  maximum  r.  This,  however,  overestimates  the  effect  since  the  flow 
nearer  the  stagnation  streamline  is  more  effective  in  transferring  momentum 
to  the  sphere.  On  the  basis  of  this  premise  and  the  fact  that  most  of 
the  data  is  taken  farther  downstream  and/or  with  smaller  spheres,  the 
radial  variation  in  flow  properties  is  neglected.  Radial  impact  surveys 
by  Ashkenas5’  also  indicated  a  negligible  variation  of  Mach  number  or 
impact  pressure  over  a  Tew  nozzle  diameters  about  the  axis  for  large 
nozzle-probe  separations. 
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The  effect  of  axial  gradients  on  the  drag  is  somewhat  more  difficult 
to  assess.  One  criterion  is  that  —  /  p  be  negligible  over  d  .  This  may 
be  too  stringent  a  requirement  since  in  hypersonic  low  density  flow  only 
the  front  half  of  the  sphere  is  effective.  In  any  event,  the  neglect  of 
axial  gradients  wren  using  the  larger  spheres  may  be  questioned.  The 
situation  may  be  further  complicated  when  considering  nearly  free  molecular 
flow.  Since  the  molecules  reflected  from  the  front  of  the  sphere  "see" 
an  increasing  density  gradient  upstream,  perhaps  more  of  Ihem  are  scattered 
back  onto  the  surface  than  would  be  the  case  in  a  uniform  stream.  On  the 
other  hand,  radial  gradients  would  tend  to  remove  more  molecules  from  the 
area.  The  results  reported  here  are  based  either  on  the  density  of  the 
sphere  center  location  or  the  stagnation  point  and  axial  gradients  are  not 
o+herwise  accounted  for, 

E.  Effects  of  Temperature  and  Background  Gas 

The  effect  of  "freezing"  of  the  temperature  and  Mach  number  at  the 
lower  values  of  P  d  on  the  experimental  drag  coefficient  is  negligible 
and  hypersonic  conditions  are  maintained.  There  are  some  problems  due 
to  the  appearance  of  +he  free  sfream  temperature  and  speed  ratio  in  the 
theoretical  expressions  and  these  effects  are  considered  when  they  occur. 

A  most  important  consideration  in  using  the  free  jet  as  a  flow  field 
is  the  effect  of  the  background  gas  and  shock  system  on  the  core  flow. 

These  effects  have  been  investigated  by  a  number  of  authors^  with 
particular  emphasis  on  criteria  for  maintaining  an  isolated  core  flow. 

It  has  been  maintained  that  the  Mach  disc  location  can  be  calculated, 
however  the  extent  of  its  influence  upstream  of  this  "position"  due  to 
thickening  at  low  background  pressures  is  not  sc  clear. *  The  oblique 
side  shocks  must  also  broaden  as  the  density  is  lowered  and  may,  at  some 

X 

‘See  for  example  References  4,  19,  20,  and  21. 

*Bier  and  Hagena19  indicated  that  the  Mach  disc  became  diffuse  enough  to 
disappear  from  flow  photographs  as  a  Knudsen  number  defined  as  Kn  =  A  /Y 
approached  unity.  c  m 


IB 


point,  influence  the  centerline  properties.  At  low  enough  background 

pressures  the  jet  may  become  "porous"  to  background  gas  and  the  core 

flow  can  be  Influenced  by  the  chamber  pressure.  Although  the  results 

of  these  effects  and  the  values  of  parameters  characterizing  them  are 

perhaps  not  satisfactorily  determined  in  general;  it  is  clear  that  these 

effects  are  to  some  degree  a  function  of  P  /P  and/or  the  absolute  value 

o  c 

of  the  background  pressure.  Since  the  cere  flow  must  be  dependent 
of  the  chamber  pressure  for  the  assumptions  and  calculations  used  to  be 
valid,  the  drag  force  must  not  be  influenced  by  a  reasonable  variation 
in  background  pressure.  The  chamber  pressure  was  varied  in  the  tests 
reported  here  and  data  affected  by  such  a  variation  will  be  noted  and 
considered  to  be  in  error. 

Finally  the  possibility  of  condensation  can  be  considered.  Molecular 
beam  experiments  indicate  that  no  condensation  occurs  at  the  conditions 
encountered  for  nitrogen.  For  argon  some  condensation  may  be  present 
but  no  effect  on  the  data  was  observed. 
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SECTION  iV 

EXPERIMENTAL  RESULTS 


A.  General  Discussion 

Drag  data  taken  at  axial  positions  from  10  to  22  nozzle  diameters 
downstream*  are  normalized  in  the  form  of  a  drag  coefficient  and  pre¬ 
sented  as  a  function  of  a  number  of  familiar  flow  variables.  The  selection 
of  the  independent  parameters  was  influenced  by  comparisons  with  avail¬ 
able  experimental  and  theoretical  results  as  well  as  an  attempt  to  determine 
a  representa+ive  variable  for  drag  in  transition  flow  that  could  pe  cal¬ 
culated  with  reasonable  accuracy.  Variation  of  Ihe  nozzle  size  and  the 
axial  position  of  the  sphere  allowed  an  increased  range  of  transition  flow 
to  be  investigated  and  provided  assurance  that  the  results  were  not 
specifically  related  to  the  nozzie  or  jet  geometry.  Three  sphere  sizes 
were  used  to  cover  a  wider  range  of  the  independent  variable.  The  flow 
properties  are  calculated  at  the  location  of  the  center  line  except  where 
otherwise  notes. 

The  results  are  divided  roughly  Into  three  groups  corresponding  to 

X  XX 

small  '  10-13),  intermediate  -  ~  13-18)  and  large  (~  >  22)  nczzle- 
n  n  n 

sphere  separations.  However,  this  division  is  not  rigorously  maintained 

and  some  overlap  will  be  obvious  in  the  presentation.  The  selection  is 

arbitran  ly  made  to  show  the  agreement  of  the  data  taken  at  these  stations. 

A  difference  was  evident  when  comparing  the  data  taken  close  to  the 

nozzle  and  that  taken  at  the  largest  distances  downstream.  At  fixed 

Knudsen  number,  the  drag  coefficient  data  is  systematically  higher  for 

stations  close  to  the  nozzle,  and  decreases  with  increasing  x  /d  .  This 

trend  is  reduced  by  calculating  the  flow  properties  at  the  position  of 

the  sphere  stagnation  point  on  the  jet  axis  instead  of  the  center  location. 

The  results  of  this  analysis  are  given  i r,  sections  V  and  VI. 


One  additional  station  was  +aken  at  54  nozzle  diameters  with  the  0.445 
mm  nozzle. 


20 


B.  Drag  Measurements  with  Source  and  Model  Initially  at  Room  Temperature 

The  drag  coefficient  for  nitrogen  is  presented  as  a  function  of  the 
Knudsen  number  (assuming  pX-const.)  for  low  to  high  vaiues  of  nozzle-  . 
sphere  separation  in  Figures  4.1  through  4.3.  Figures  4. A  and  4.5  pre¬ 
sent  similar  results  for  argon.  For  aii  data  the  source  temperature  was 
about  298°K.  Theoretical  inviscid  Mach  numbers  ranged  from  8.7  to  I?. 8 
for  nitrogen  and  from  15.0  to  25.6  for  argon. 

Data  outside  of  the  vertical  lines  shown  on  the  figures  are  no+ed  to 
be  in  error  and  the  effect  of  the  background  pressure  variation  on  these 
data  is  given.  Although  these  limits  could  often  be  derermined  by  examin¬ 
ing  the  graphs,  they  were  obtained  experimentally  in  each  case  by  recording 
any  change  in  force  with  independent  variation  of  the  chamber  pressure  at 
constant  source  conditions.  The  position  of  the  vertical  line  excludes 
the  values  of  Knudsen  number  at  which  the  smallest  recordable  change  in 
force  (<  0.3$  of  sphere  weight)  occurred  with  a  reasonable  increase  in 
background  pressure  (obviously  an  unlimited  increase  in  always  produced 
an  effect).  Although  such  a  determination  is  admittedly  qualitative,  the 
agreement  with  the  inflections  observed  in  the  drag  curves  indicates  that 
the  test  is  valid.  The  limiis  were  recorded  for  each  sphere  size  and 
occurred  at  approximately  the  same  source  pressure  for  a  given  nozzle- 
sphere  separation.  Clearly  these  limits  appear  at  different  Knudsen 
numbers  as  a  function  of  sphere  site  but  only  those  for  maximum  Knudsen 
number  with  the  0.635  mm  sphere  and  minimum  Knudsen  number  with  1/8  inch 
spheres  are  shown.  All  intermediate  points  affected  are  omitted  for 
clarity  since  a  valid  point  could  be  obtained  at  that  Kn  under  another  set 
of  conditions.  One  adidtional  verification  of  the  ‘'unaffected'*  data  is 
provided  by  The  fact  that  the  curves  for  different  sphere  sizes  are  in 
agreement  and  in  many  cases  an  area  of  overlap  is  shown.  The  dal  a  in¬ 
fluenced  by  the  background  pressure  variation  do  not  agree  in  this  manner. 

The  exact  nature  of  the  influence  of  jet  shock  structure  and  back¬ 
ground  pressure  on  the  core  flow  has  not  been  quantitatively  resolved. 
Measurements  of  sphere  drag,  however,  do  not  represent  a  satisfactory 
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method  for  determining  the  details  of  these  effects.  The  sensitivity  of 
the  balance  does  provide  an  accurate  indication  of  the  onset  of  such 
effects  on  the  drag  measurements  and  the  fact  that  no  sting  extended 
through  the  shock  structure  assured  that  the  effect  was  not  precipitated 
by  the  influence  of  a  mechanical  support. 

At  high  densities  the  sphere  bow  shock  may  be  well  formed.  This 
shock,  however,  becomes  a  broadened  density  disturbance  as  the  density 
is  reduced.  The  jet  Mach  disc  and  oblique  side  shocks  also  broaden  as 
the  pressure  is  lowered.^  Interactions  between  these  shock  structures, 
the  sphere  wake,  and  background  gas  appear  to  account  for  the  effects 
observed.  When  the  Knudsen  number  is  high,  the  influence  of  background 
pressure  causes  an  increase  in  drag.  It  is  difficult  to  see  how  this  can 
be  attributed  to  the  broadening  Mach  disc  since  one  would  expect  a  de- 
crease  in  impact  pressure  as  the  tail  of  the  shock  is  entered.  Perhaps 
an  'nfi  Itration  of  gas  through  the  side  shocks  and  accelerated  by  the 
outer  part  of  the  jet  increases  the  centerline  impact  pressure  at  lew 
densities.  Another  possibility  would  be  tne  influence  of  the  continually 
broadening  oblique  side  shocks  contributing  to  This  increase.  These 
effects  could  not  be  distinguished  by  the  techniques  of  this  experiment. 

The  influence  of  background  pressure  on  the  measurements  at  high 
source  pressures  was  a  decrease  in  drag.  Gregorek  and  Luce20  have  shown 
that  the  pressure  on  the  downstream  side  cf  the  sphere  can  be  increased 
with  little  effect  on  the  upstream  face  as  the  data  is  taken  further  down¬ 
stream  (closer  to  the  Mach  disc).  The  effect  occurs  even  though  the  ncmina 
position  of  the  Mach  disc  is  a  number  of  nozzle  diameters  downstream  of 

'See  references  5,  19,  20. 

f 

In  fact  this  effect  was  observed  during  the  tests.  For  a  number  of  runs 
the  measurements  are  taken  at  continually  reduced  source  pressure  and 
correspondingly  reduced  background  pressure  due  to  the  characteristics 
of  the  system.  The  drag  coefficient  increased  smoothly  to  a  maximum 
and  then  decreased.  Credit  for  the  decrease  is  given  to  entering  the 
Mach  disc  shock  structure. 
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tne  sphere  position.  This  appears  to  De  borne  out  by  the  decrease  in 
drag  observed  at  lew  Knudsen  numbers  here.  The  decrease  occurs  later 
(at  both  lower  Knudsen  numbers  and  higher  source  pressures)  if  the 
sphere  is  positioned  further  upstream.  Grer  'ek20  shows  that  the  sphere 
bow  shock  can  actually  be  deflected  by  the  je,  shock  structuie. 

It  is  not  implied  ihat  the  detailed  nature  of  the  interaction 
mechanisms  are  proved  by  the  resulls  presented  here,  but  only  tnat  these 
effects  could  account  for  the  results  observed.  More  detailed  analysis 
and  careful  experiments  would  be  necessary  in  an  attempt  to  make  use  of 
the  data  influenced  by  the  jet  structure  and  no  such  effort  has  yet  been 
made.  The  data  affected  is  noted  to  be  in  error  and  shown  both  to  indicate 
that  care  was  taken  to  assure  the  validity  of  the  results  presented  and 
that  such  an  investigation  is  essentia!  when  using  the  free  jet  as  a 
flow  field. 

The  free  molecular  limits  shown  in  Figures  4.1  -  4.?  are  calculated 
f  rorrrz 


00 


with 

c  =  o'  -  !  diffuse  reflection 

o  =  o'  =0  specular  reflection  . 


Since  the  stagnation  temperature,  the  original  sphere  temperature, 
and  the  surrounding  room  temperature  were  the  same  for  these  tests,  the 
sphere  is  assumed  to  be  at  stagnation  temperature  for  the  calculation. 

T 

An  exception  will  be  noted  as  T  insulated. 

w 
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The  data  generally  indicate  that  this  assumption  is  valid.  The  diffuse 
reflection  limit  was  assumed  since  all  of  the  data  rise  above  the  specular 
limit  and  indicate  that  the  diffuse  value  is  approached.  Molecular  beam 
experiments  with  a  target  of  the  same  material  as  the  sphere  models  also 
indicates  diffuse  reflection  at  room  temperature.'  Although  the  theore¬ 
tical  values  of  the  speed  ratio  and  free  stream  temperature  used  may  be 
in  error  due  to  freezing  at  large  nozzle-sphere  separations  and  low  source 
pressures,  calculations  based  on  roughly  estimated  frozen  values1*4'23  in¬ 
dicated  that  the  effect  was  negligible  for  the  calculation  of  Cpp^. 

The  most  attractive  feature  of  Kn  as  a  flow  parameter  for  corre¬ 
lating  the  data  ir  that  it  can  be  accurately  calculated.  Values  for 
viscosity  (and  hence  mean  free  path)  are  not  available  at  the  local  free 
stream  temperatures  encountered.  However,  since  the  molecules  reflected 
from  the  sphere  surface  have  energies  correspondi ng  to  stagnation  tempera¬ 
ture,  it  appears  that  the  effect  of  molecular  collisons  causing  a  departure 
from  free  molecular  conditions  might  be  betier  character i zed  by  a  Knudsen 
number  based  on  source  conditions  than  a  parameter  dependent  on  free 
stream  properties.  The  data  are  presented  as  a  function  of  other  parameters 
elsewhere  in  this  chapter. 

A  small  amount  of  scatter  in  the  results  .  hown  in  Figures  4.1  -  4.5 
is  evident.  The  data  includes  all  points  for  each  run  except  where  points 
are  obscured  by  the  density  of  the  plot.  No  averaging  or  selection  was 
employed  to  reduce  the  scatter.  The  apparent  scatter  is  largely  due  to 
the  presentation  of  data  at  different  nczz ie-sphere  separations  on  the 
same  plot.  The  runs  were  repeated  and  reproducibility  was  better  than  the 

4. 

scatter  shown  on  the  graphs. L  No  effort  was  made  to  label  separate  runs 
at  the  same  conditions,  and  often  only  one  run  is  shown. 

^Taken  by  the  Gas-Surface  Studies  Laboratory,  Department  of  Aerospace 
^Engineering  and  Engineering  Physics,  University  of  Virginia. 
xThe  one  exception  is  the  data  at  12  mm  where  some  repeatability  problems 
were  encountered  due  to  nozzle  vibration.  However,  the  run  is  representa¬ 
tive  and  the  difference  was  at  worst  -  5%  with  runs  taken  over  the  span 
of  a  year. 


29 


C.  Drag  Measurements  with  Source  and  Model  Heated 

A  limited  amount  of  data  was  taken  with  a  heated  source  using  the 
0.95  mm  nozzle  and  the  results  are  shown  in  Figure  4.6.  The  purpose  of 
source  heating  was  to  increase  the  local  temperatures  so  that  free 
stream  properties  could  be  calculated  and  with  more  confidence  to  move 
the  Mach  disc  downstream.  However,  no  significant  gain  was  achieved 
at  the  source  temperatures  attained. 

Due  to  experimental  problems  only  the  1/8  inch  spheres  were  used. 

The  sphere  temperature  is  now  known  and  the  nozzle  position  may  have 
changed  with  increasing  temperature.  Some  decrease  in  the  drag  coefficient 
might  be  expected  if  the  sphere  was  cooled  by  radiation  to  the  ambient 
surroundings  but  no  such  decrease  was  observed.  Radiation  from  the  heated 
nozzle  may  have  compensated  somewhat  for  this  cooling.  If  the  sphere  was 
near  the  nozzle  stagnation  temperature,  no  effect  of  source  heating  on  the 
drag  coefficient  would  be  expected  unless  the  surface  accommodation  model 
changed.  The  parameter  Kn  may  not  be  valid  -  correlating  this  data  if 
the  surface  is  cooled,  however,  the  results  remain  slightly  higher  than 
the  ambient  data  in  terms  of  Kn^  or  Re2-  Finally,  the  sphere  magnetic 
moment  and  hence  the  calibration  could  be  affected.  It  was  concluded 
that  no  measurable  change  in  drag  coefficient  due  to  source  heating  could 
be  determined  from  these  results. 

One  positive  result  is  indicated  by  the  curve.  The  maximum  rise  in 

the  data  due  to  jet  influence  is  reduced  by  heating  the  source.  This 

I  /4 

might  be  expected  since  the  flow  field  size  increases  as  (T  )  ,  indi¬ 

cating  that  heating  the  source  should  provide  an  increased  data  range. 
However,  better  control  over  the  heating  process,  considerably  higher 
temperatures,  and  careful  calibration  would  be  essential  for  any  signifi¬ 
cant  improvement. 

Attempts  were  made  to  heat  the  sphere  surface  and  monitor  the  surface 
temperature  using  infrared  techniques.  The  results  are  shown  in  Figure 
4.7.  Again  the  variable  Kn  may  not  be  the  correct  correlating  parameter 
for  these  data  if  the  scattering  is  diffuse.  However,  qualitatively  the 


30 


31 


FIGURE  4.6  Source  Temperature  Effect  on  Sphere  Drag 
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FIGURE  4.7  Drag  Coefficient  for  Heated  Spheres 


increase  observed  indicates  diffuse  reflection  since  no  change  would 
be  expected  if  the  scattering  were  specular.  These  data  would  be  more 
valuable  near  free  molecular  flow  but  such  results  were  not  successfully 
obtained.  Calculation  of  the  increase  in  drag  expected  in  free  molecular 
flow  for  completely  diffuse  reflection  indicates  that,  at  least  when 
correlated  in  terms  of  Kn,  the  measured  increase  in  due  to  surface 
heating  is  approximately  the  same  or  greater  than  the  calculated  increase 
(AC^  '  0.2  @  2I5°C).  This  measured  difference  would  be  less  if  some 
increased  effective  Knudsen  number  were  used  that  included  the  effect 
of  the  increased  body  temperature. 


0.  Discussion  of  Experimental  Error 
An  estimate  of  the  error  in  the 
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measured  quantities  is  as  follows: 

+  0.05  torr 
±0.5  torr 
±  0.01  mm 
±  0.04  mm 
±  0.005  mm 
±  0.5°C 
±  0.0013  mm 
±  0 . 00 i 3  mm 
±  0.0025  mm 
±  0.05  mm 
±  \% 

±  0.25? 

±  0.15? 

-  ±  0.25?  of  sphen  support 

current 

<  *  0.2 % 

-  ±  I  % 

±  0.5% 
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The  error  resulting  from  the  use  of  the  theoretical  equations5  is: 

P./P  ±  \% 

i  o 

M  ±  \/2% 

If  the  above  errors  are  all  added,  the  maximum  possible  error  is 
~  4. 5-6. Q%  (from  low  to  high  Knud sen  numbers)  in  the  drag  coefficient 
at  the  smallest  nozzle  sphere  separation.  The  error  is  less  for  larger 
distances  downstream.  Typical  maximum  error  flags  are  shown  in  Figure 
4.1  and  4.3.  Repeatability  of  the  data  is  much  better  thsn  these 
maximum  possible  error  values. 

The  error  due  +o  axial  and  radial  gradients  in  the  flow  field  has 
been  discussed.  The  radial  gradients  are  considered  negligible  and 
the  axial  gradient  problem  is  believed  to  be  a  minimum  when  the  flow 
properties  are  calculated  at  the  sphere  stagnation  point  location. 

The  error  due  to  gradients  is  considerably  less  for  large  nozzle-sphere 
separations,  indicating  a  preference  for  these  data. 
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COMPARISON  WITH  RESULTS  OF  OTHER  EXPERIMENTS 

A.  Flow  Properties  Based  on  Sphere  Center  line 

The  current  literature  on  the  measurement  of  sphere  drag  at  low  densities 
contains  dala  obtained  by  a  number  of  different  techniques  and  representing 
a  wide  range  c i  key  variables  such  as  the  Mach  No.  and  Knudsen  No.  In  order 
to  make  a  meaningful  comparison  of  these  data  some  correlation  parameter 
must  be  used.  Although  the  literature  also  contains  much  discussion  con¬ 
cerning  the  proper  choice  of  correlation  parameter,  one  that  seems  to  have 
received  relatively  wide-sprea  use  is  the  Reynolds  No.  downstream  from  a 
normal  shock  as  computed  by  isentropic  theory.  We  have  arbitrarily  chosen 
this  as  the  independent  variable  or  correlation  parameter  and  have  denoted 
i t  as  Re^. 

In  Figure  5.1  the  results  of  the  previous  section  are  presented  as  a 
function  of  Re^,  where  the  dependent  variable  is  the  drag  coefficient 
normalized  to  the  free  molecular  diffuse  value  with  the  sphere  at  stagnation 
temperature.  This  presentation  is  also  consistent  with  that  of  many  other 
authors.  Also  shown  on  this  figure  are  the  free  flight  data  of  Kinslow 
and  Potter24  for  nitrogen,  Bailey25  for  air  and  the  free  flight  shock- 
tunnel  results  of  Geiger26  for  air.  All  of  the  results  of  Kinslow  and  Potter 
a^e  presented.  Bai iey's  data  continues  to  much  higher  Reynolds  numbers  than 
are  within  the  range  of  the  graph.  The  reduction  in  scatter  with  the 
magnetic  suspension  technique  is  evident  and  can  be  evfen  less  if  data  at 
onlv  one  sphere  station  is  presented.  An  extension  of  the  data  toward  free 
molecular  flow  is  provided  by  the  present  technique  and  although  it  appears 
that  the  free  molecule  limit  is  attained,  the  data  to  the  left  of  the  vertical 
line  are  of  questionable  validity  due  to  the  influence  of  the  jet  structure. 
This  again  points  out  the  need  for  caution  when  using  the  jet  flow  field. 

All  of  the  free  flight  results  are  for  relatively  cold  wall  conditions. 
This  suggests  that  the  present  results  should  perhaps  be  greater  than  the 
free  flight  data.  This  increase  is  not  apparent  in  the  comparison.  However, 
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FIGURE  5.1  Normalized  Drag  Compared  with  Other  Experimental  Results 
(Flow  Properties  at  Sphere  Center  Location) 


The  present  data  definitely  tend  toward  the  free  molecular  limit  for  stagna¬ 
tion  wall  temperatures  in  all  tests. 

In  section  IV-A  mention  was  made  of  the  fact  that  the  drag  coefficient 
tended  to  be  higher  when  measured  at  smaller  nozzle-sphere  separations. 

The  maximum  difference  observed  was  about  8-10$  when  comparing  data  at 
constant  Knudsen  number  for  1/8  inch  spheres  at  the  minimum  and  maximum 
X  /d  values  for  nitrogen.  The  difference  dec. eases  with  increasing 
Knudsen  number  to  about  4$  at  near  free  molecular  conditions.  Any  such 
discrepancy  was  considerably  less  when  comparing  any  other  stations  and 
always  less  for  argon. 

B.  Flow  Properties  Based  on  Sphere  Stagnation  Point 

Although  there  is  some  slight  Mach  number  effen  (the  free  molecular 
drag  coefficient  is  greater  at  lower  Mach  numbers),  it  does  not  appear  to 
account  for  the  difference  observed.  Even  for  the  smallest  distances  at 
which  tests  were  made  the  flow  was  hypersonic  and  very  little  Mach  number 
effect  is  expected  in  this  regime.  A  more  reasonable  assumption  is  that 
this  difference  is  due  to  the  axial  gradients  in  the  flow  which  cause  a 
greater  error  for  stations  close  to  the  nozzle  and  for  larger  spheres. 

The  difference  trend  was  eliminated  for  the  1/8  inch  spheres  by  calculating 
the  flow  properties  at  the  location  of  the  stagnation  point  of  the  sphere 
instead  of  the  center  line  location. 

Some  slight  difference  as  a  function  of  nozzle-sphere  separation  re¬ 
mained  at  higher  Knudsen  numbers  (smaller  sohere  models)  and  this  was 
attributed  to  the  mass  flow  corrections  discussed  in  C  apter  II.  The 
mass  flow  corrections  were  greater  at  the  lower  source  pressures  used  to 
obtain  high  Knudsen  numbers  and  even  greater  at  the  further  reduced  source 
pressures  to  produce  the  same  Knudsen  number  at  smaller  nozz le-sDhere  separa 
tions.  Therefore,  if  the  data  is  slightly  over-corrected  by  the  effective 
nozzle  size  used,  this  effect  would  result  in  the  difference  observed.  In 
any  event,  any  consistent  rise  in  the  drag  coefficient  remaining  with  de¬ 
creasing  ^s/dn  is  always  less  than  3$  for  data  not  influenced  by  the  jet 
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flow  field.  This  difference  must  be  considered  to  be  within  the  error  of 

a 

the  results  in  view  of  the  assumptions  used  in  the  data  reduction. 

Data  reduced  using  the  flow  properties  at  the  location  of  the  sphere 
stagnation  point  are  shown  in  Figure  5,2,  Although  th's  improved  the 
internal  agreement  between  the  data  taken  at  different  sphere  stations, 
the  comparison  with  the  free  flight  results  is  less  impressive  at  higher 
Reynolds  numbers.  The  disagreement  is  confined  to  the  1/8  inch  spheres 
and,  as  has  been  mentioned,  these  are  the  least  credible  of  the  present 
data.  Since  the  sphere  size  is  approximately  3  times  the  nozzle  diameter, 
the  axial  gradient  influence  is  not  negligible.  It  would  seem,  however, 
that  the  stagnation  point  is  the  more  representative  location  from  momentum 
transfer  considerations. 

The  presentation  of  data  in  terms  of  properties  at  the  front  of  the 
sphere  indicates  the  same  smooth  increase  toward  the  free  molecular  limit 
shown  previously.  These  results  are  given  in  Figure  5.3  for  nitrogen  and 
Figure  5.4  (low  X^/d^)  and  Figure  5.5  (high  /dn )  for  argon  as  a  function 
of  free  stream  Knudsen  number. 

An  additional  comparison  is  made  with  the  results  of  Smolderen,  Wendt, 
and  Naveau.27  These  results  were  taken  with  a  sting-mounted  sphere  in  a 
free  jet  using  air  at  Mach  8.3  (X/dn  =  9.1).  The  comparison  is  shown  in 
Figure  5.6.  There  is  good  agreement  with  the  present  nitrogen  data  taken 
at  low  nozzle-sphere  separations  and  based  on  sphere  center  line  properties 
However,  the  data  taken  further  downstream  and/or  based  on  stagnation  point 
properties  predict  a  later  rise  to  the  free  molecular  limit  than  Smolderen' 
results.  (This  can  be  seen  in  Figure  6.9  of  the  next  section.)  There  may 
be  some  Mach  number  effect  when  comparing  the  data  with  those  at  large 
nozzle-sphere  separations  arid  the  possibility  of  sting  effects  cannot  be 
excluded.  Also  the  calculation  of  Kn^  is  questionable  for  the  present 
resu I ts. 
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FIGURE  5.4  Sphere  Drag  vs.  Knudsen  No.  (Low  xs/dn)  for  Argon 
(Flow  Properties  at  Stagnation  Point  Location) 
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FIGURE  E.6  Sphere  Drag  vs.  Knudsen  No.  (Low  xs/dn)  for  Nitrogen 
(Flow  Properties  at  Sphere  Center  Line) 


Sphere  drag  data  has  been  taken  by  other  investigators  at  low  Mach 
numbers'1"  in  transition  flow  but  those  presented  inciude  the  available 
hypersonic  results.  No  comparison  was  available  for  argon. 


*1* 

See  for  example  references  28  and  29. 


SECTION  VI 

COMPARISON  WITH  THEORETICAL  MODELS 


Comparison  will  be  made  with  the  results  obtained  by  Willis30*31 
using  the  modified  Krook  model,  the  linearized  Krook  solution  of  Rose,32 
and  the  first  collision  results  of  Baker  and  Charwat.33  Few  theoretical 
computations  for  sphere  drag  In  near  free  molecular  flows  are  available 
and  these  are  believed  to  represent  those  somewhat  applicable  to  the 
present  results. 


A.  Near  Free  Molecule  Approaches 

Willis31  obtained  a  solution  for  sphere  drag  using  the  modified  Krook 
model  of  Boltzmann's  equation,  assuming  the  body  to  be  moving  in  a  gas  in 
thermal  equilibrium  at  T^  with  no  ablation  or  absorption  and  the  molecules 
leaving  the  body  to  have  a  Maxwellian  di stri butiv.i  corresponding  to  some 
temperature,  T^.  The  analysis  predicts  a  reduction  of  the  free  molecular 
drag  in  terms  of  the  small  parameter 

d  I  ,L 

a  =  6n»r*(? vh  (6J) 

b 

where  l/Sn^  is  a  relaxation  time  for  translational  energy  and  T^  is  the 
temperature  characterizing  the  reflected  molecules.  The  constant  6  was 
determined  from 

kT 

6  -  — —  .  (6.2) 


If  it  is  then  assumed  that  the  viscosity  can  be  obtained  from  the  Chapman- 
Enskog34  relation 


=  PAS 

irm 


1/2 


(6.3) 


for  hard  sphere  molecules,  the  parameter  a  can  be  written  as 


r  T  1/2 

a  =  (11)  C,^)  /  Kn„. 


(6.4) 
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Using  these  relations  an  analytic  approximation  to  the  computer  solution 
for  the  drag  of  a  sphere  was  given  by  Maslach,  et  a  I.33  as 


f-'or  the  calculation,  equation  (6.5)  was  rearranged  in  the  form 


(6.8) 


T  T 
00  00 

Using  theoretical  values  of  *— ■ 1  =  ■*»—  and  S,  =  1.871  the  resul  rs  for 

To  'b  b 

nitrogen  are  shown  in  Figure  6.1  and  6.2.  For  the  high  ^s/dn  results 
shown  in  Figure  6.3  considerable  liberty  is  taken  in  estimating  the 
frozen  value  of  from  Marrone14  for  the  theoretical  curves  (to  predict 
a  maximum  reduction  from  free  molecular  flow  with  the  theory)  and  better 
agreement  was  obtained.  Marrone's  results  are  for  rotational  freezing 
and  are  higher  than  the  freezing  values  for  translation  indicated  by 
molecular  beam  studies.36  Argon  results  are  shown  in  Figure  6.4  for 
high  X  /dn  and  =  1,581.  Figure  6.5  for  argon  at  low  X^/d^  essentially 
repeats  the  data  at  high  nozzle-sphere  separations.  The  agreement  of  the 
data  with  the  modified  Krook  model  is  good  even  though  the  condition 


FIGURE  6.1  Sphere  Drag  vs.  Free  Stream  Knudsen 
Nitrogen 
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FIGURE  6."  Sphere  Drag  vs.  Free  Stream  Knudsen  No.  (Int.  x  /d  )  for 
Nitrogen  s  n 
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FIGURE  6.5  Sphere  Drag  vs.  Free  Stream  Knudsen  No.  (Low  x$/dn)  For 
Argon 


>  2  is  net  fulfilled  for  the  test  results.  The  agreement  !s  surprising 
when  the  assumptions  made  to  calculate  Kn^  are  considered.  The  solution  is 
only  valid  for  small  a  and  hence  will  fail  at  low  Knudsen  numbers.  This 
is  evident  in  Figure  6.2. 

Baker  and  Charwat33  obtained  a  first  correction  to  the  free  molecular 
f'ow  over  a  sphere  by  assuming  that  only  first  collisions  between  molecules 
are  importanl.  It  was  further  postulated  that  t-he  sphere  was  in  a  parallel 
and  uniform  beam  of  billiard-ball  molecules,  that  reflected  molecules  have 
an  average  speed  much  less  than  free  stream  molecules,  and  that  the  mean  free 
path  of  reflected  molecules  is  independent  of  molecular  velocity.  The  as¬ 
sumption  of  "slow"  reflected  molecules  is  violated  in  the  present  tests  and 
would  appear  better  for  cold  wall  conditions.  Willis31  rearranged  the  resulms 
of  Baker  and  Charwat  in  terns  of  the  parameters  used  here  and  the  resulting 
equation  is 

Cn  =  cr._,  -  (0.24S,  +  1.06)  /  Kn  .  (6.9) 

D  DFM  b  ® 

The  results  were  shown  in  Figures  6.1,  6.3  and  6.4  and  are  seen  to  be  in 
serious  disagreement  with  the  data. 

Rose32  used  the  original  Krcok  equation  with  the  sphere  represented  by 
the  addition  o*  a  point  source  term.  The  equation  was  iinearized  about  the 
distribution  function  at  infinity  and  solved  using  Fourier  transform 
techniques.  The  assumptions  used  for  the  point  source  model  and  the 
linearization  are  diffuse  reflection  and 

d  <<  D  <<  X 
s 

and  D  is  the  observer  position  where  the  solution  is  valid.  This  restricts 
the  solution  to  near  free  molecular  flow.  The  first  order  solution  obtained 
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(6. 10) 


CD  =  CDFM 


2 

Kn, 


1/2 


0.33  (y-0 


-  0.17  /(y) 


1/2 


M 


where  KnD  is  a  parameter  used  in  the  model.  Rewriting  equation  (6.10)  in 
R  t 


terms  of  the  variable  used  hero 


CD  '  CDFM 


(fi  j |  (0.33  - 


0.  12, 


(6.11) 


The  results  are  ^hown  in  Figures  6.1  through  6.4,  The  data  are  seen  to 
depart  quickly  from  agreement  with  these  results. 

Sherman  et  a  I.37  postulated  that  for  free  jet  testing  the  comparison 
with  theory  might  best  be  made  in  terms  of  a  parameter  based  on  source 
conditions.  This  is  indicated  since  such  a  oarameter  can  be  readily 
calculated  and,  with  the  body  at  about  stagnation  temperature,  the  im¬ 
portant  irttermolecular  collisions  may  be  characterized  by  this  temperature. 
The  parameter  suggested  is 


p  V  d  T  1/2 
00  °°  s  ,  o.  ,  „c 
a  =  —  \v—)  /  4S 

o  V  T  /  O 

o  b 


(6.12) 


where 


mV  2 


S  2  =  — 
o  2kT 


(6.13) 


Here  we  will  assume  that  T^  =  Tq  and  the  parameter  cxq  is  related  to  a  given 
in  equation  (5.6)  by 


T  v 
(— )ct 
T~  yo 


(6. 14) 


Equations  (6.8),  (6.9)  and  (6.10)  can  be  rewritten  in  terms  of  a  and, 


Maslach  et  al.j5  listed  these  values.  Slightly  different  constants  result, 
depending  on  the  relationship  used  to  relate  KnR  and  Kin,*.  However,  the 
comparison  with  the  data  is  qualitatively  the  same  in  all  cases. 
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interpreting  a  as  a 

Cn  =  CnrM  -  a 
D  DFM  o 

CD  =  CDFM  ‘  ao 
LD  =  CDFM  ‘  “o 


the  resulting  equati 


0.37  t 


3.32  _  2.57 

Sb  S.  2 

d  J 


0.75  - 


0.26 


I  + 


4,36" 

Sb  J 


ons  are 

(Wi I  I  is) 


(Rose ) 


(Baker  and  Charwat) 


(6.15) 


(6J6) 


(6. 17) 


Data  reduced  using  properties  at  the  front  of  the  sphere  are  compared 
with  equations  6.15  -  6.17  in  Figures  6.6  and  6.7  for  nitrogen  and  Figure 
6.8  +or  argon.  Here  we  see  a  clear  preference  for  agreement  with  the 
modified  Krook  solution. +  The  droop  occurring  at  low  values  of  is  due 
to  data  influenced  by  the  jet  structure.  With  nitrogen  the  sphere  was 
assumed  to  be  at  stagnation  temperature  for  the  calculation  of  Cpp^  and 
agreement  is  good.  The  argon  data  are  in  better  agreement  if  the  equilib¬ 
rium  temperature  for  an  insulated  sphere  is  used  to  calculate  the  free 
molecular  value.  In  any  event,  the  choice  of  Cpp^  does  not  change  the 
slope  of  the -curve.  The  data  shown  in  Figure  6.7  are  taken  further  down¬ 
stream  and  the  boundary  layer  correction  is  less.  This  would  appear  to  be 
the  best  data  and  good  agreement  is  indicated  with  the  linear  dependence 
pred i cted  by  Willis.31 


The  data  can  be  compared  again  in  terms  of  the  free  stream  variable 

Kn^.  This  comparison  is  shown  in  Figures  6.9  (low  X^/d  )  and  6.10  (high 

X  /d  ).  Although  the  agreement  is  slightly  worse  in  terms  of  stagnation 
s  n 

point  p>  nerties,  the  modified  Krook  Model  is  clearly  indicated  to  be  the 

most  applicable.  Due  to  the  fact  that  a  can  be  calculated  more  accurate  I v 

o 

and  would  seem  to  better  represent  the  collison  process,  preference  is  given 
to  the  linear  comparison. 

Data  based  on  center  line  properties  show  the  same  agreement  with  the 
tendency  tc  be  slightly  beiow  the  prediction  of  Willis. 
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FIGURE  6.6  Comparison  with  Theory  for  Nitrogen  (Stagnation  Point 
Properties,  Int.  x  /d  )  Based  on  Source  Parameter 


1 


FIGURE  6.10  Comparison  with  Theory  for  Nitrogen  (Stagnation  Point 
Properties,  High  x$/dn)  Based  on  Knudsen  No. 


B.  Approaches  from  Continuum  or  Slip  Flow 

Davis  and  Flugge-Lo+z38  obtained  a  solution  for  departure  from  continuum 
flow  by  solving  the  boundary  layer  equations,  including  second  order  effects, 
for  hypersonic  flow  over  a  blunt  body.  Only  their  numerical  solution  ob¬ 
tained  for  a  sphere  at  Mach  10  with  a  wall  to  stagnation  point  ratio  of  0.6 

can  be  compared  here.  The  resuits  were  only  compufed  to  first  order  in 

I  /2 

inverse  (Re2)  and  are  shown  in  Figures  6.11  and  6.12.  The  solution  is 
for  y  =  1.4  and  viscosity  proportional  to  the  square  root  of  +he  temperature. 
The  value  0.89  is  used  for  the  pressure  drag  coefficient  at  high  Reynolds 
number.  Theoretical  wall  to  free  stream  temperature  ratios  for  the  data 
are  -  20  in  Figure  6.11  and  ~  25  in  Figure  6.12.  These  compare  to  an  assumed 
value  of  12.6  for  the  model.  The  model  must  fail  at  low  Reynolds  numbers 
and  is  only  shown  to  indicate  that  an  extension  of  the  data  would  be  in 
general  agreement  as  illustrated  by  Figure  6.12.  The  disagreement  in  Figure 
6. II  can  be  removed  by  basing  the  flow  properties  on  the  stagnation  point 
values,  thereby  lowering  the  high  Reynolds  number  portion  of  the  curve.  The 
effect  of  this  correction  would  be  less  in  Figure  6.12.  Figure  6.13  gives 
the  argon  data  in  terms  of  the  familiar  parameter  of  Reynolds  number  behind 
a  normal  shock.  The  behavior  is  similar  to  that  for  nitrogen. 

Aroesty28  suggested  the  same  form  for  representing  the  drag  coefficient 
in  transition  flow,  i.e.. 


Kinslow  and  Potter28  determined  that  K2  was  negative  by  fitting  their 
data  fo  this  equation.  The  inviscid  pressure  drag  term  has  been  determined 
by  Hodges39  to  be  0.92  and  later  by  Bailey25  to  be  0.885.  The  value  of  0.89 
(in  agreement  with  Davis38)  was  assumed  here  for  nitrogen  and  the  dafa  for  a 

Using  the  freezing  value  of  rotational  temperature  from  Marrone14  and 
average  values  of  P  d^  for  the  estimate,  the  ratios  would  be  reduced  to 
~  II  and  ~  15  resDecti ve  ly . 
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FIGURE  6.12  Sphere  Drag  vs.  Reynolds  No.  Behind  a  Normal  Shock 
(Int.  x  /d  )  for  Nitrogen 
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FIGURE  6.13  Sphere  Drag  vs.  Reynolds  No.  Behind  a  Normal  Shock. 
(High  *s/dn)  for  Argon 


number  of  runs  were  f i tied  by  lessl  squares  to  equation  (6.18).  Typical 
results  of  the  fittings  are  given  in  Table  2.  In  all  cases  K2  is  negative. 
The  range  of  the  values  is  evident  as  a  function  of  nozzle  sphere  separation. 
Note  that  these  constants  are  overly  influenced  fcy  a  large  amount  of  data 
at  very  low  Reynolds  numbers  and  not  sufficiently  balanced  by  high  Reynolds 
number  data.  The  estimated  value  of  0.88  was  used  for  with  argon. 

The  values  spf  Kj  for  nitrogen  can  be  compared  with  the  values  of  Davis38 
of  3.8  for  ~  12.6  and  2.8  for  ~  =  4.2  at  Mach  10. 


TABLE  2 

CONSTANTS  FOR  DATA  FIT 


Nitrog  n 

x  /d 
s  n 

K1 

k2 

Standard 

Deviation 

o 

10.9 

3.9 

-2.  1 

8 . 7 ( I0~2) 

13.6 

3.1 

-1.5 

5. i ( 1 0~2) 

15.8 

3.0 

-1.4 

1 „5( I0~2) 

16.4 

2.9 

-1.3 

3.  1 ( I0"2) 

Argon 

13.6 

3.6 

-1.7 

3,8(1 0~2) 

16.4 

3.  1 

-1  .1 

3.6( I0"2) 

21.8 

3.  1 

-1.3 

2.5( I0‘2) 
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SECTION  VI  I 


SUMMARY  AND  CONCLUSIONS 

c/'lensive  drag  measurements  have  been  made  on  spheres  In  the  transition 
regime  of  gas  dynamics.  The  data  exhibit  reduced  scatter  and  improved 
repeatability  compared  with  results  previously  available.  The  techniques 
employed  have  provided  sting-free  measurements  while  maintaining  control 
over  the  model  in  a  high  velocity  flow.  A  wide  range  of  Knudsen  numbers  in 
near-free  molecular  flow  has  been  investigated  for  both  nitrogen  and  argon 
gases. 

The  findings  of  the  investigation  are: 

1.  The  combination  of  the  magnetic  balance  and  free  jet  techniques 
represents  a  productive  and  comparatively  inexpensive  method  for 
obtaining  accurate  drag  data  In  high  speed  low-density  flow. 

2.  General  agreement  has  been  demonstrated  with  the  relatively 
limited  data  available  for  comparison. 

3.  The  measurements  verify  the  expected  smooth  increase  in  drag  through 
the  transition  regime  tending  asymptotically  to  the  free  molecular 
limit  (based  on  a  diffuse  thermally  accommodated  reflection  model). 
There  is  no  true  tendency  for  any  of  the  data  to  rise  above  this 

I imit. 

4.  Good  agreement  Is  obtained  over  a  reasonably  wide  range  of  flow 
conditions  with  the  theoretical  results  of  Willis  using  a  modified 
Krook  model . 

5.  Studies  made  with  temperature  variation  of  both  the  source  gas  and 
the  sphere  surface  offer  qualitative  agreement  with  expectations; 
however,  the  experimental  technique  has  not  been  sufficiently  re¬ 
fined  to  allow  quantitative  comparison. 

In  addition,  several  observations  can  be  made  concerning  the  use  of  the 
free  jet  as  a  low  density  flow  field. 
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1.  Correlation  parameters  in  terms  of  the  source  conditions  and  the 
local  density  are  those  most  accurately  determined  in  free  jet 
+esting  and  they  appear  to  provide  a  successful  correlation  of 
the  variation  of  drag  coefficient  in  transition  flow. 

2.  Drag  measurements  at  large  nozzle-sphere  separations  are  the  most 
repeatable,  exhibit  the  least  scatter,  and  require  the  least  cor¬ 
rection  for  effective  nozzle  size  at  the  seme  local  flow  conditions. 
The  latter  is  important  since  accurate  mass  flow  measurements  (used 
to  correct  the  nozzle  size)  are  difficult  to  obtain.  However, 
attention  must  be  paid  to  additional  aspects  such  as  temperature 
and  Mach  number  "freezing"  when  testing  at  large  nozzle-sphere 
separations. 

3.  The  influence  of  the  jet  shock  structure  and  background  gas  on  the 
sphere  drag  may  be  successfully  monitored  by  the  chamber  pressure 
variation  tests. 

4.  For  a  fixed  Dumping  system,  low  nozzle-sphere  separations  are  best 
for  testing  at  the  near-continuum  end  of  the  transition  regime, 
while  near  free  molecular  conditions  are  better  simulated  at  large 
nozzle-sphere  separations. 

5.  The  discrepancy  between  data  obtained  at  small  nozz ie-sphere 
separations  and  those  obtained  at  large  separations  (occurring 
primarily  for  larger  spheres  in  nitrogen)  can  be  reduced  by  cal¬ 
culating  the  local  flow  properties  at  the  leading  stagnation  point 
rather  than  at  the  sphere  center. 
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